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1 . SUMMARY 


A theoretical study was performed on the effects of 
gravity level and environmental composition upon the struc- 
ture of laminar hydrocarbon diffusion flames. The study 
was made on vertically oriented, axisymmetric, fully developed 
laminar fuel jets issuing into a quiescent atmosphere. 

The mathematical formulation includes the effects of 
inertia, viscosity, diffusion and gravity. Combustion is 
treated with an equilibrium model, and both Fick's Law and 
multicomponent representations were used to model the 
diffusion process. These mechanisms were coupled in a bound- 
ary layer type formulation and the resulting equations were 
solved by an explicit finite difference technique. 

The principal data used for comparison with the model 
predictions was obtained from a NASA Lewis drop tower. Flame 
shapes for methane under normal gravity and zero gravity 
conditions were provided by color photography. The Reynolds 
number based on fuel jet conditions ranged from approximately 
5 to 300. The remaining data were obtained from the literature 
and were limited to normal gravity conditions. A dimensional 
analysis shows that the flame structure is characterized by 
the maximum flame width radius, velocity and thermodynamic 
state. With the aid of the theoretical model these proper- 
ties have been related to their counterparts at the jet exit 
plane . 

The comparison between theory and experiment for steady 
state flames showed excellent agreement for all normal gravity 
flames while under zero gravity conditions the predictions 
at the lower Reynolds numbers showed generally wider flames 
than were observed. The absence of kinetics and radiation 
effects was shown to be the primary" cause for this difference 
between theory and experiment. The addition of partial oxi- 
dation models in the formulation substantiated this conclusion 
while the experimentally observed "orange-reddish" colors lent 
credence to the presence of soot and "cool" or slow oxidation, 
zones within the flames . 



2. INTRODUCTION 


The structure of laminar diffusion flames has -been the 
subject of investigation for a great many years. Applica- 
tions include en ergy utilization with interest expanding 
more recently to include the general problem of fire hazards. 

It is the latter problem that we are most interested in here • 

In this connection the ignition and propagation of the fires 
will invariably occur in a quiescent or very low velocity 
environment where the basic combustion process does indeed . • •••-. , . 
involve laminar flow. Moreover, our specific interest is in- 
fires aboard spacecraft where the effects of gravity level 
and environmental composition in terms of oxygen concentra- .... 
tion are crucial additional considerations. ,, 

A review of the existing literature shows that in each 
case a significant gap exists between the particular set of 
experimental data and the isolated correlations that have been 
attempted. The gap is in establishing a basic understanding, 
of the structure of laminar diffusion flames in terms , of first 
principles. That is, there has been a dearth of systematic • 
mathematical modeling on this problem. 

The recent work of References 1 and 2 has adequately y. 
summarized the state-of-the-art regarding such mathematical 
modeling of laminar diffusion flames. The existing work may 
be categorized in terms of (1) semi-empirical approaches based 
upon correlation of specific data typified by the work of 
References 3 and 4; and (2) special solutions of the conserva- 
tion equations typified by References 5 and 6. 

The semi-empirical approaches were ground-breaking works 
but nevertheless produced correlations which differ in func- - 
tional dependence between the various authors . Thus, while 
attempting to provide an understanding of the problem, a degree 
of confusion and misleading conclusions have resulted. This 
is typified by the recent work of Reference 7, where a combina- 
tion of attempts to extend existing correlations with the aid ; 
of euristic arguments concludes that the laminar diffusion ' \-u 
flame does not depend explicitly upon a Reynolds number . The 
result of this is a prediction of flame length as a function of 
gravity level which is just the opposite of what is in fact 
observed ! 
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Common to the special solutions of the conservation 
equations is the absence of the gravity force while the 
assumption of similar solutions is also made. Specifically, 

Burke and Schumann, Reference 5, developed a classical solu- 
tion for a ducted diffusion flame assuming uniform constant 
velocity and constant coefficients. This reduced the problem 
to one of solving the single homogeneous linear partial 
differential equation for an element mass fraction character- 
izing the fuel. The solution of such an equation is relatively 
simple but it does not contain the mechanisms that are 
necessary to characterize the true structure of diffusion 
flames. More importantly, the benefit of obtaining a simple 
solution in this manner precludes systematic extensions 
necessary to build a more realistic model. Reference 6 
treats the unbounded laminar diffusion flame and solutions 
are obtained assuming similarity while gravity effects are 
neglected. Again, the limitations are analogous to these 
associated with the Burke -Schumann problem. 

What has been required is an approach which at the start 
recognizes the mechanisms which are of potential importance 
while providing a basis necessary to readily and systematically 
extend the model to include additional mechanisms. Specifically, 
the coupling of velocity, thermal and chemical species fields, 
is needed. Furthermore, the mechanisms which influence their 
behavior must be coupled together within an analytical framework 
of sufficient generality so that such primitive assumptions 
including similarity and constant transport properties are not 
dictated by the solution technique and as such are not required. 

This report is devoted to the development of a mathematical 
model and solution technique, which is designed to provide the 
generality needed to gain a basic understanding of laminar 
diffusion flames influenced by gravity. The approach involves 
a finite difference solution of the describing partial differen- 
tial equations for mass, momentum and energy applied to a basic 
vertical axisymmetric fuel jet configuration. 

This approach recognizes that the velocity, temperature and 
concentration fields are dependent upon each other while the 
finite difference solution eliminates restrictions on the gener- 
ality of the boundary conditions as well as providing a frame- 
work which may be readily extended without requiring' a new' solu- 
tion technique to be developed from scratch. The effects * included 
are: diffusion, inertia, viscosity , combustion and gravity. 
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The work described in this report includes the mathematical 
formulation , a comparison of the theory with experiments, a 
dimensional analysis , and a parametric study- In addition, a 
critical evaluation of the weaknesses of the model is given and 
a delineation of model refinements is presented. 


3„ ANALYTICAL MODEL FORMULATION 


A schematic of the flow field of the vertical laminar jet 
is shown in Figure 1- The flow is axisymmetric and for now is 
viewed in a steady state configuration. The fuel jet velocity 
profile will generally be nonuniform and for fully developed 
laminar flow will be parabolic. The pure fuel case involves a 
flat ' (constant) concentration profile across the jet.: Similar- 

ly, for an adiabatic tube (T Q = T e ) the fuel temperature will be 
constant across the tube. Thus, even if the mechanisms forcing 
the respective development of the velocity, temperature and 
species, fields were similar, there is non-similarity in these 
boundary conditions ruling out, in general, similar solutions. 

Now, the velocity, temperature, and species fields are 
controlled by the distribution of momentum, energy and mass. 

Momentum . - The fuel is injected at some finite velocity 
into the quiescent surroundings. The velocity difference pro- 
duces a shear force resulting in an interchange of momentum 
between the fuel jet and its surroundings. This tends to equal- 
ize the velocities which in our case means a deceleration of 
the fuel velocity toward zero. This process of diffusive trans- 
port of momentum causes the jet to spread and entrain environ- 
mental gas into the mixing region. Simultaneously, a pressure 
nonuniformity develops due, in general, to the velocity varia- 
tions and due to any environmental pressure variation impressed 
upon the mixing zone. This may be "favorable" (causing accelera 
tion) or "adverse" (causing deceleration). In the jet problem 
it is instructive to consider the "impressed" pressure variation 
only, assuming that the static pressure is radially invariant 
and the flow is boundary-layer like. In this case the local 
pressure gradient (in the axial direction) is equal to the 
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pressure gradient in the environment* This is a favorable press-* 
ure gradient since it is equal to the gravity force (-gp) • The 
combination of pressure gradient and local gravity force yields 
the ’’buoyant" force and will cause an acceleration of the flow 
for (p ~/o) > 0. Note, however, that the value of the buoyant 
force can be locally negative* Thus, in the near jet fuel rich, 
cool region, away from the flame front, the "buoyant" force 
will be negative for heavier-than-air fuels such as butane* 

This can result in stagnating the flow locally with the attendant 
formation of a recirculation zone. This will be discussed later 
in connection with a comparison of our model with Wohl 1 s experi- 
ments (Reference 4) . For light fuels like methane, the buoyant 
force is always positive. This is of particular interest in 
connection with the recent work, previously mentioned, by Cochran 
and Masica, References 1 and 2. 

Thus, there are three direct effects upon the velocity fields 
(a) shear or molecular diffusion of momentum, (b) pressure grad- 
ient, and (c) gravity, where (b) and (c) combine to give the 
buoyant force in the jet problem* 

Energy . - The temperature field is affected by several pro- 
cesses. In the vertical jet problem the convection of energy 
is a primary mode of energy transport and therefore depends upon 
those mechanisms which alter the velocity field. In general the 
velocity field can alter the thermal field directly through vis- 
cous dissipation and by the rate at which non-dissipative pressure 
work is done. These effects are negligible in our low speed 
problem. Additional mechanisms in the jet problem which will be 
shown to be most important are the heat release rate due to com- 
bustion, the heat conduction away from the flame front, and the 
energy transported due to diffusion of species. Thus, the 
temperature field also depends upon composition in terms of its 
effect upon density, specific heat, the transport coefficients as 
well as the direct effect upon the species production rate, 
which is an exponential function of temperature. 

Species Diffusion . - The concentration field is affected by 
the convection of species through the flow and therefore by the 
velocity field. The mechanisms which directly "force" changes 
in a given specie are the chemical production rate of that specie 
and the diffusion of that specie* As the fuel comes into contact 
with the environmental gas the resulting concentration gradients 
produce a diffusional flux of environmental gas into the fuel (and 
fuel diffuses outward into the environment) . These diffusional 
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gradients produce an air/fuel ratio distribution ranging from 
fuel rich on the axis (x £ L, Figure 1) to fuel lean in outer 
regions of the mixing zone. The surface upon which the equiv- 
alence ratio is unity will correspond to the maximum tempera- 
ture in the flow field and is referred to as the "flame front." 

The primary assumption is i that diffusive processes are 
important only in the radial direction. 

Governing Equations ,. - The development of an analytical 
tool can serve in several ways: (1) provide direct predictions 

(2) interpret experimental data and (3) provide the basis for 
meaningful dimensionless parameters which characterize the 
process and as such can aid in organizing experimental data. 

The first two aspects are an inherent part of a formulation 
written directly in terms of the variables of the problem in 
dimensional form. To aid in achieving the latter aspect from 
a basic formulation it is helpful to introduce appropriate 
characteristic; quantities :and thereby non-dimensional ize the 
describing equations. It is customary practice to do this by 
normalizing the variables directly in terms of the boundary con 
ditions. However, it will be shown here that the resulting 
dimensionless groups do not truly characterize the process and 
that this is due to the large variation in flow properties 
throughout the domain of interest. Accordingly, the current 
approach involves the introduction of a set of characteristic, 
quantities which are based upon some state within the flow that 
more accurately reflects the flame structure. Although this 
study has shown that the state of flow at the maximum flame 
width point (including the flame radius and velocity) is appro- 
priate, the development given here is general. For this 
purpose the following variables are introduced: 


u s u/u* 
r s r/r* 
p * p/p* 
H = H/H* 
V = v/u* 


x = x/r* 

U = 


A0 


p s -p 



p e - p 

p b -p* 


(i) 


Where, as cited above, the starred quantities are evaluated 
at a point in the flow which best characterizes the process. 
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Assuming (1) that diffusive processes are important 
only in the direction normal to the primary flow direction, 
(2) that combustion is diffusion controlled, and (3) that 
the flow is steady and axisymmetric , the describing equa- 
tions in non-dimensional form are: 

Continuity 


d (pur) 5 ( pvr ) 

dx dr 


0 


( 2 ) 


Momentum 


pu I# + 9V f§ = (^) i i tiiai/ari 
dx dr Re r dr 


(3) 


+ (-~) Ap 

Res' 


Diffusion of Elements 




— bay — ba 
pu T"=r + p\r ~r= 
dx dr 


_ *1., 1 

v Re y r dr 


(4) 


Energy 


— dH — dH 
pu 7" + pv — 
dx dr 


, 1 x _1 d r yx_ dh 

Re r dr i Pr dr 


- r £ (h k j k )| (5) 
k J 


where the above parameters: Reynolds, Grashof numbers are 

defined ash ■ • . 


Re = 


p*u*r* 

ji* 


Reynolds No. 


Gr _ iilEi 

12*2 P* 


= Grashof No. 


( 6 ) 


To complete the description of the problem in terms of 
the basic variables , ( u , T and g^) , additional relations are 

required. Specifically, the equations of state, expressions 
for the diffusional fluxes, j and j , and the definition of 
the chemical state in terms of a combustion model, are needed. 
These are described below: 

Equations of State 

P = — z — ZT7 ~T ‘ i = all gaseous species (7) 

RT £ OC'/ w 1 

i 

h = £ aV (8) 

k 

where 

h = H - u^/2 (9) 


Diffusional Fluxes, j and j 

For laminar flows both Fick 1 s Law and multicomponent models 
have been employed to express the diffusional fluxes in terms of 
the basic properties of the flow. The Fickian model is simple 
and can be expected to be adequate under conditions where the 
diffusing gases do not differ significantly in their respective 
molecular weights. The multicomponent representation is alge- 
braically complex and expresses the direct coupling of each 
diffusional flux to the entire species field. This model is 
potentially more accurate than Fick's Law, depending upon how 
well the set of binary diffusion coefficients are known. 


a. Fickian Diffusion 

Assuming that all diffusion coefficients, A are equal, 
is equivalent to reducing the mixture to a set of J binary sub- 
systems comprised of the specie of interest and the remaining 
mixture. The diffusional mass flux for the k*-* 1 specie is then 
given by: 
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( 10 ) 



For this case it is possible to define a Schmidt number such 
that 

J Sc dr 


or 


.k u* - da 

3 = Scr* M Sr 


and 


■rk 1 - da 

- 1 Sc ^ ar 


( 12 ) 


The corresponding diffusional mass flux for each £ element 
is given by: 


r r\ t 


1 

1 - da 


3 Sc ** Sr 


(13) 


where 


2- a k 
k w* 


(14) 


so that 


~l - w .k 

3 “ f * "T 3 

k w 


(15) 


and v is the amount of element- l in specie k. Now, substi- 
tution of Eqs . (12) and (13) into the above conservation equations 

introduces the Schmidt Number as an additional parameter, where 
Sc = ratio of diffusivity of momentum to the diffusivity of mass. 
Thus, the working forms of the element and energy conservation 
equations are given by: 
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( 16 ) 


and 


'•'-t 

t- da . 9a 

Pu — + pv 
dx dr 


ReSc 


ii_,- do\ 

r dr dr 


— dH — dH lid 
Pu — - + pv ~ ^ = — -• “ 
dr dr Re r or 


;{ 


i“tfc 


dh 

dr + 


Sc 


__ -k 
S h 


da 

dr 


Si 


(17) 


As part of the "simple" transport mechanism model the 
Schmidt Number and Prandtl Number are treated as specified 
quantities. This reduces the problem of transport property 
definition to that of specifying a relationship between the 
absolute viscosity, /i, and the local thermodynamic properties. 
For this purpose the Sutherland Law is used based upon the 
properties for air: 




X .458 x 10 6 T 1 ' 5 . kg 

T + 110.4 m-sec. 


(18) 


where T is in °K. Although this expression is for air 
(independent of the actual composition) the primary dependence 
upon temperature is retained. The differences that do exist 
between the levels of viscosity predicted by the Sutherland 
Law and the actual values associated with, say, pure fuel, air 
and products of combustion are found to have some effect upon 
the accuracy of the flame structure predictions but not to 
the extent of altering the conclusions regarding the mechanisms 
responsible for the particular flame behavior. Nevertheless, 
to provide a quantitative measure of the potential improvement 
in the accuracy of the predictions, a detailed transport pro- 
perty model including multicomponent diffusion, has been 
developed. 
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b. Multicomponent Diffusion 


As cited above, potential improvements in the accuracy of 
the predictions are afforded by a detailed transport property 
representation. The improvements can be expected to be most 
significant at low Reynolds Number and at low g-levels where 
diffusive processes dominate, and for chemical systems where 
the dominant species vary significantly in their respective 
molecular weights. 


For non-equal diffusion coefficients the generalized multi- 
component diffusion formulation given in Reference 8 is appro- 
priate. In terms of mass fractions, cfc, the required relations 

are given by: 


N 

- E 4 — 

p j-i 


/ k . j 

(a j J 


j .k X 

) 


N 

£ 


1 , da 

— (a . 

. - w. u 
3=1 3 


k 


dr 


- a 




dr 


(19) 


and 


N 

r j 3 = o 

j = 1 


( 20 ) 


where N is the number of species. Equation Q-9) provides N-l 
independent relations for the j diffusional mass fluxes and 
Eq. (20) completes their description. The diffusional mass 
fluxes for the elements are related to corresponding species 
fluxes, those relations of the type given by Eq. (15), viz. 



k£ w 


w 


.k 

J 


( 21 ) 


The details of the matrix inversion technique and the 
coupling to the conservation equations are described in 
Appendix A. To accomplish the desired solution certain infor- 
mation is required by Eqs . (19) and (20). Specifically, the 

mass fractions, their gradients and the density must be known. 
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These are, in fact, determined in the course of generating 
the flow field as discussed in detail in Appendix A. In 
addition, the must be specified and these are determined 

based upon the kinetic theory analysis given in Reference 8 , 

viz ; ■ ' : •' 1 


= 2.663 x 10 ~ 2 T 1 ’ 5 . ,m , 

Kj P O (1,1) r kj sec 

kj 

o N 

Where [T] = K, and [p] = 

In addition, 

i / \ +w i ' 

" V J 2 \ w j 

with the collision cross-section given by: 


( 22 ) 


(23) 


a + a . 
_k 1 


kj 


(24) 


Finally, the collision integral, Q , y . , as tabulated in Refer 
ence 8 as a function of the reduced ^temperature, T/€^., has 
been curve fit for the present application. The resulting 
equation is given by; 




(1?D = i . 5146 - 0.62499 T . + 0.10023 t. 2 . (25) 


kj 


kj 


Where, 


and 


T. . = -tn (T/c . ) 
k] k3 


(26) 


«/. = <A,e . 

k] k 3 


(°K) 


(27) 
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The required properties for the participating species are 
given in Table I, below. 


Now, consistent with this detailed computation of the 
diffusion coefficients, both the viscosity, fi , and the con- 
ductivity, X, are computed using the following mixture rules. 
Reference 8. 


Viscosity 

N gases 

jLZ = £ 

k-1 




■ [ 1+ ^ ft ] S 

CL + Z ■ i -“ — ^a. 

. , y. . 3 

3=1 *3 

jA 


(28) 


where 


^3 


w. h 
fc 1 ) 




k 


(29) 


and 


7 *3 


w. 


w. 


k 7 


w. 

8 (1+ — ) 
v W. 

D 


(30) 


with 




2.6693 x 10~ 6 /tW* 


k a k 8 \(2,2) 


(JSO_) 

m-sec 


(31) 


and. 


= 1.6507 - 0.6688 T k + 0.10725 T^. 2 


(32) 


(2 0 ) . 
where Q. v ' has been curve fit from the tabulation given 

in Reference 8. 
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TABLE I 


CONSTANTS 

FOR TRANSPORT 

PROPERTY 

CALCULATIONS 

GAS 

W 

<7 (A) 

e (°K 

H 2° 

18.016 

2.641 

809.1 

H 2 

2.016 

2 .827 

59.7 

°2 

32. 

3.467 

106.7 

CO 

28.011 

3.69 

91.7 

C0 2 

44.011 

3.941 

195.2 

n 2 

28.016 

3.798 

71.4 

CH * 
4 

16.043 

3.758 

148.6 

C (gas) 

12.011 

3.385 

30.6 


* Constants for other fuels are available for input. 
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Thermal Conductivity 


by; 


The mixture rules for the thermal conductivity are given 


X = 


where 


N gases 
v 

k=l 


\K. 


N 

a v + E a- 3 n. . C . 

k j-1 


+ 


a, X" 

k k 


N in. 

id 




j=i 

j^k 


Joule 


(m~sec-°k) 


(33) 


with 


K = 3 - 75 t^ 

k 


(34) 


and 


Jr* 

[0.352 C W“ - 0.88] 

>k = ^ f 

R 


(35) 


K = 


8.3143 x'lO 3 .: — ^ , °k 
kgm-mole 


(36) 


and 


1 + 


•fin . = 


K_ 

, X' . 

V 1 


P 




kj 


(37) 


W 


C kj * 1 + 


2.41 (W. -W.) (W, -0.142 .) 

3S 1_ k ll 

s 

(w, + w . ) " 

k 3 


Note that 


~k _ k l W 
k W 


(38) 
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Co Turbulent Transport 


Although the emphasis in the current study is on the 
structure of laminar flame s , certain of the available experi- 
mental data extends into the transitional and fully developed 
turbulent regimes. References 4 and 9« In order to provide 
a means of interpreting some of these data, provision has 
been made in the model to account for turbulent exchange of 
mass, momentum, and energy 0 The approach here involves replac- 
ing the laminar transport properties with their turbulent 
counterparts o In this regard the Fickian framework is employed 
and with the specification of ““turbulent" Schmidt and Prandtl 
numbers (usually taken equal to unity) the problem reduces to 
modeling the eddy viscosity,, Based upon previous work. 

Reference 10, the following model has been found to be adequate § 




kr^(pu) 

•3 £ 


(39) 


where k is an empirical constant and ri is a measure of the 

o "2 

scale of turbulence and is taken at the radius where the mass 
flux, (pu ) , is at its average value » 

Applications of the Pick’s Law, multicomponent diffusion 
and turbulent property transport representations are discussed 
in the subsequent sections in connection with the predictions 
of flame structure,, 


d„ Combustion Mod,el 


An assumption which has been employed in simple treatments 
of diffusion flames is that all the heat is released on the 
surface where the fuel/air ratio is stoichiometric (Reference 6) „ 
This "flame sheet" assumption is extremely restrictive and is 
not necessary. Its use is generally coupled to similar solu- 
tion techniques where accuracy is sacrificed strictly for the 
sake of mathematical simplification 0 A more accurate combus- 
tion model is offered in terms of basic equilibrium theory » 

In a diffusion flame the application of this theory leads to 
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a "local" equilibrium state at each point within the flow 
field. Thus, a solution of the conservation equations, 

Eqs „ (2), (3), (4) and (5), provides the velocity, element 
mass fractions and total enthalpy at a known local pressure 
level. Of this information the element mass fractions, 
static enthalpy and pressure, define the chemical equilibrium 
state. This includes the mass fractions and static tempera- 
ture. Because the flow field is nonuniform the local point- 
by -point equilibrium condition is referred to as "shifting" 
equilibrium. The advantages of using shifting equilibrium 
are that, (1) it is not complex, and (2) dissociation can 
be accounted for. This technique has already been success- 
fully applied in connection with application to diffusion 
flames (Reference ICL) . 

Now, in considering the potential adequacy of an equil- 
ibrium combustion model, the basic consideration was the 
residence time within the diffusion flames. It is shown 
later that the characteristic residence times for the data 
analyzed in this report are typically of the order of 10 
milliseconds . 

Experience at GASL with hydrocarbon-air finite rate 
chemistry models (c.f. References 1JL and 12) verified that 
the above characteristic times are generally in excess of 
the combustion times required for hydrocarbons at the temper- 
atures and pressures encountered in the flame. For example. 
Figure 2 of Reference 11 shows the ignition delay and total 
reaction times for a methane -air combustion process where it 
is important to note that the ignition delay and reaction 
times for most gaseous hydrocarbons, and certainly for hydro- 
gen, are shorter than that of methane. Hence, it would 
appear to be adequate to employ an equilibrium chemistry 
model to describe most normal gravity flames and stable zero 
gravity flames of interest in this report. Of course, the 
logic leading to this conclusion is based upon (1) that the 
temperature levels are, in fact, of the order of 2000 K, and 
(2) that the reaction times for "hot" combustion in premixed 
systems are meaningful for diffusive flames. 
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The basic model employed in this work was developed 
from the results of detailed equilibrium computations for 
hydrocarbon combustion. The features of the model are 
described in Appendix B„ The adequacy of the model in re- 
presenting detailed equilibrium predictions is summarized 
in Figures 3 and 4. For our current applications , typified 
by the lower initial temperature. Figure 3, the agreement is 
quite good over the entire equivalence ratio range. Although 
not crucial to the current study, it should be noted that for 
higher initial temperatures, dissociation effects become 
significant and the agreement is relatively poor in the near 
stoichiometric regime. Figure 4. 


d. Boundary Conditions 


The following conditions represent the boundary 
conditions required to set the problem for the steady axi- 
symmetric vertical jet flow in an infinite environment: 


0 £ r £ r . 

J 


x — 0 


H = H (r) 
o 

u = u (r) 
o 

a. = a. (r) 

i io 

v = 0 


Notes Arbitrary initial 
profiles may be specified 
including parabolic for 
the velocity. 

= 1 for pure fuel 


(40) 
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fu =0 for the quiescent environment 
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(41) 
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(43) 


(44) 


and the h. 's are determined with specification of the 
temperature. In addition to the boundary conditions, the 
"impressed" pressure field must be specified. This is given 
by “ 



p e g 


(45) 


In summary, the model accounts for the cross-coupling 
between the velocity, temperature and species fields as 
influenced by gravity; arbitrary transport coefficients which 
depend upon the local state; and shifting equilibrium combust io n. 


4. APPLICATIONS OF THE MODEL 


a. Fickian Diffusion 

The preceding discussion concentrated on the elements of 
a rather detailed model for diffusion flames. There are virtually 
no restrictions on the specification of boundary conditions and 
a general treatment for the transport coefficients of mass, 
momentum and energy is retained. 
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This section describes the utility of the model to 
(1) predict experimental observations, (2) delineate the 
relative importance of those mechanisms which are included 
in the model, (3) provide the basis for establishing mean- 
ingful characteristic parameters to aid in data interpre- 
tation and correlation, and (4) pinpoint theoretical and 
experimental aspects which need further consideration* The 
bulk of our attention during this study has been devoted to 
the experimental work for which both normal gravity and zero 
gravity data was available - namely the methane-air flames 
of Cochran and Masica (c.f. References 1, 2 and data supplied 
in a private communication from Mr. Thomas H. Cochran, 

June 23, 1971). To supplement this work, the normal gravity 
data of Wohl, et al. (Reference 4), and Hawthorne and Hottel 
(References 3 and 9) was examined. Also, parametric studies 
were made for methane-air, ethylene-air and propylene-air 
flames, and the results are available for comparison with 
future experimental work. 

The first step in the application of the model was to 
verify that we could successfully predict the normal gravity 
flame lengths of Reference 1. As shown in Figure 5, good 
agreement was obtained between the analysis and the data for 
these relatively large fuel nozzle radius flames by postul- 
ating a Prandtl number of 0.8 and a Schmidt number of 0.5. 

It may be noted that when unity Prandtl and Schmidt numbers 
were assumed (for the initial conditions of Test 20, Refer- 
ence 1) , a theoretical flame length of twice the experimental 
length was obtained. When the higher fuel flow rate/smaller 
nozzle experimental data became available, Xprivate communica 
tion from Mr. Thomas H. Cochran, June 23, 1971), a further 
comparison of theoretical and experimental flame lengths was 
made, as shown in Figure 6. It was found that for the higher 
fuel flow rates (Q > 5 cm^/sec) that the use of a Schmidt 
number value of 0.6 yielded better agreement with the experi- 
mental data. This value was used in all subsequent calcula- 
tions. For the initial zero-g- comparisons the operating 
conditions corresponding to the preceding normal-g methane 
flames were used. A comparison of flame length predictions 
with the experimental results for a series of normal-g and 
zero-g • ; methane cases is summarized in Figure 7. 
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These direct comparisons show generally excellent agree- 
ment for the normal-g flames and fair agreement for the steady 
state 0-g flames. In particular , for these specific sets of 
methane cases the predicted 0-g flame lengths are generally 
somewhat shorter than the observed lengths. In all cases, how- 
ever, the observed trends for this set of data are predicted 
where here the 0-g flames are longer than the norma 1-g flames. 
In addition, these buoyant flames are too complex to be 
treated by similar solutions, and thus analyses such as that 
by Goldburg and Cheng (Reference 6) are inappropriate. For 
example, the non-similarity of the calculated radial velocity 
profiles for a methane-air normal gravtiy flame is shown in 
Figure 8. In the near, cool inner region where diffusion 
dominates, the flow is decelerated. Then, downstream as the 
flow is heated due to combustion, buoyancy becomes important 
and the flow is accelerated. 

Nevertheless, the preceding study for these 0-g flames 
showed that the predicted flames, based upon Fickian diffu- 
sion, were not only shorter but broader than indicated by 
the experiments, particularly at the lower Reynolds numbers 
as shown on Figure 9. 


b. Multicomponent Diffusion 

The theoretical calculations described in the preceding 
discussions were based upon the Fickian representation for 
diffusive transport. As cited in Section 3, the importance 
of detailed transport property considerations suggested the 
inclusion of a multicomponent diffusion model and comparisons 
have been made with the Fickian diffusion model and with 
experimental data where it was available. These comparisons 
are shown in Figures 10 through 13. In general it appears 
that for the hydrocarbon flames there are not gross differ- 
ences between the multicomponent and the Fickian analyses. 
Figures 10 through 12 . However, for the hydrogen flame 
there does appear to be a significant effect. Figure 13. 

This is to be expected since the diffusion coefficients are 
sensitive to the molecular weights of the chemical species. 
However, comparison was made using an isolated piece of normal 
gravity laminar data of Hawthorne (Reference 9) , and shows 
that the Fickian flame length (0 = 1) is about 150 cm while 
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the multicomponent prediction yields a length of between lib 
and 120 cm. Although the multicomponent prediction is closer 
to the experimental value, L ^ , the agreement is still poor. 
Although this prediction shows the potential importance of 
multicomponent diffusion, the discrepancy between the pre- 
diction and the single data point needs further experimental 
information before this limited comparison can be truly 
assessed. 

Furthermore, the potential importance of soot in the 0-g 
hydrocarbon flames dictates that a multicomponent representa- 
tion be used to account for slow diffusion of the soot, 
particularly relative to the gaseous species. 

In general, then, it is desirable to retain the multicom- 
ponent transport model for further comparisons and general 
refinement of the flame model as well as for further studies 
in connection with sooty flames. 


c. Turbulent Diffusion 


Although turbulent flames are not of direct concern in 
this study, certain available data have been analyzed to 
demonstrate the overall utility of the model for analyzing 
turbulent as well as laminar flames. For this purpose the 
turbulent hydrogen and propane normal gravity flame data of 
Hawthorne and Hottel (Reference 3) was used. The turbulent 
viscosity model used, as discussed above, was of the form: 

m t = k r yp u) g 

As shown in Figures 14 and 15, good agreement with experimental 
results was obtained for both the hydrogen and the propane 
flames. 
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Thus, it appears that the same model framework may be 
used for laminar and turbulent flames provided that the 
appropriate transport property relations are employed. 


d. Negative Buoyancy and Plow Reversal 


In terms of overall data evaluation an interesting and 
rather important observation was made. In particular, a 
start was made in analyzing the butane-air data of Wohl, et 
al. (Reference 4) as reported on pp. 58-60 of Reference 13. 

It was found there that the experimental flame length for 
the highest fuel flow rates (Re Q =3,140) could be repro uce 
by imposing a viscosity level, typical of transitions 
(laminar-to-turbulent) while for a lower fuel rate (Re Q z 
the flow field calculation was terminated by the occurrence 
of a stagnation point indicating the start of a recirculation 
region. The recirculation is due to the persistence , for the 
heavier than air fuel, of a negative buoyancy effect m the 
near jet region which is equivalent to a strong adverse 
pressure gradient* 


Upon examining the butane data of Wohl in more ^detail , 
it was found that recirculation was predicted even for fuex 
flow rates as high as 100 cm 3 /sec, with initially turbulent 
fuel Reynolds numbers (Re 0 = 2,030)* Velocity profiles 
illustrating the onset of the recirculation region are shown 
in Figure 16. This is a very clear indication that the 
possibility of local recirculation for heavier than air fuels 
must be kept in mind, and anticipated for other "heavy hydro 

carbons . 


As will be discussed below, negative Grashof number pro- 
files were also obtained in the near region of the propylene 
evaluation case (where the fuel/air weight ratio was 1.59), 
although acceleration due to heating occurred prxor to the 
onset of recirculation* The fact that the analysis of u3.1 ° 
the laminar butane-air data was sensitive to and significantly 
affected by the negative buoyancy phenomena, suggests additional 
areas of importance in terms of practical applications* T is 
includes fires not only involving fuel vapors which are heavier 
than air but also conditions where the flame is decelerated* 
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Thus, further exploration of the occurrence of recirculation 
is indicated. 

The preceding comparisons between theory and experiment 
show excellent agreement between the basic Fick's Law model 
and the normal gravity data, while for low Reynolds number, 
zero gravity flames, the model yields wider steady state 
flames than observed. However, the inclusion of a multi- 
component model for diffusion provided some reduction in 
this discrepancy. Finally, application of the basic model 
to turbulent flames shows excellent agreement between theory 
and experiment. 

Now, the excellent agreement obtained for the normal 
gravity flames suggested that a more detailed examination 
be made to determine the steps required to improve the model 
in the low Reynolds number, zero gravity regime. Potential 
mechanisms responsible for the discrepancy between the basic 
theory and experiment at these conditions include kinetic, 
elliptic and transient effects. 



MODEL REFINEMENTS 


5 . 


a - Kinetics Effects - Partial Oxidation 


The fact that many of the zero gravity flames were ob- 
served to have ant orange-reddish appearance, seemed to indicate 
that significant amounts of solid carbon were being formed 
by pyrolysis in the fuel rich portions of these f lames • In 
addition, such colors are typical of cool flame phenomena 
involving partial oxidation of the fuel. Both pyrolysis/soot 
formation and the cool flame oxidation process will reduce 
the temperature levels in the actual flames. Finally, the 
presence of soot in these flames will tend to further reduce 
the temperature levels through the mechanism of thermal radia- 
tion o 


One of the principal effects of a reduction in tempera- 
ture level is a decrease in the diffusion rate with an attend- 
ant reduction in flame width. To -provide a means of assessing 
the importance of such phenomena a set of simple ’’partial" 
oxidation models have been employed. Specifically, the 
partial oxidation models employed were: 

(a) Fuel is burned to H^O and C(s) 

(b) Fuel is burned to CO and (and C(s) in 
the most fuel rich portion of the flame) 

(c) Fuel is burned to: H^O and CO 

These models are "psuedo" complete combustion models selected 
such that the specification of the "element" concentration 
defines the species mass fractions without the need for any 
detailed "equilibrium" computations. As such, they are similar 
to the basic complete combustion model detailed in Appendix B. 
Of these, the first two models are the most significant since 
the first provides for the maximum production of solid carbon 
possible, and thus enables us to observe the extreme case of 
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total fuel pyrolysis without subsequent downstream carbon oxida- 
tion. The second model also provides for a significant reduction 
in the flame temperature level. Also, the limit of non-reacting, 
or frozen chemical flow, with the pure fuel mixing with the 
ambient air, was also employed to provide a "base line" (minimum 
mixing rate) for comparison purposes. 

As may be seen in Figure 17, a dramatic reduction in 
theoretical flame width was obtained by using the partial oxi- 
dation models. Indeed, the theoretical results for r m /r Q now 
lie within the scatter of the zero gravity steady-state flame 
data. More detailed information is shown in Figures 18 through 
21 . 

The methane-air flames shown in Figures 18 and 19 (tests 
7 and 8 of Private communication from Mr. Thomas H. Cochran, 

June 23, 1971), mark the transition between an extinguished zero 
gravity flame (Q = .75 cmVsec) and an established, steady zero- 
g flame (Q = 1.08 cmVsec) for .051 cm radius burner. 

Figure 19 shows the experimental data together with analy- 
t ic predictions made using three combustion models as designated 
on the figure. In addition, the stoichiometric surface is shown 
for a calculation made assuming no combustion. Here it is 
important to note that the respective maximum temperature 
"flame" surfaces shown correspond to stoichiometric conditions 
associated with the particular combustion model. The frozen 
flow stoichiometric surface, however, corresponds to the 
"complete," or full chemistry model. Now, considering the 
stoichiometric surface for the full chemistry model with and 
without combustion, it is seen that this surface becomes long 
and narrow, as the heat release is reduced. However, the stoi- 
chiometric surfaces, corresponding to the partial oxidation 
models show smaller flame radii, as desired, but also shorter flame 
lengths. Figure 19. This behavior is observed even for the 
steady state predictions of flames that actually quenched during 
the drop test sequence. Figure 18. These results suggest that 
the heat release mechanism is not a uniform one throughout the 
flow field. Rather, in the near jiet region partial combustion 
producing soot and partially oxidized species dominates the 
reaction process while upon flowing downstream the continued 
he&ting and increased residence time provides the conditions 
for these species to begin to burn off. As a consequence, the 
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actual flame front would traverse the various "fixed" stoichio- 
metric surfaces associated with each partial oxidation model 
tand ; is consistent with the experimental observations as shown, 
for example, in Figure 18. 


The extent to which our findings for the methane-air 
system can be generalized to other fuel/oxidizer/diluent systems 
needs to be determined. Ideally, the availability of zero 
gravity flame data for both hydrogen-air and high carbon to 
hydrogen ratio fuel systems, such as acetylene-air or propylene- 
air, would make it possible to isolate the contribution of the 
soot kinetics to the zero gravity flame quenching process. 

While quantitative data for the propylene-air system was not 
available during the course of this study, some qualitative 
results did become available for ethylene-air and propylene-air 
diffusion flames (c.f. Communication from Mr. Thomas H. Cochran 
to R. B. Edelman), and this encouraged us to perform certain 
parametric studies for these systems. In particular, it is 
known that unlike the methane-air flames the steady-state zero 
gravity propylene-air flames are smaller than the equivalent 
normal gravity flames. Also, it has been remarked that the 
ethylene and propylene flames (with a hydrogen/carbon molar 
ratio 6f 2) appear to be redder and sootier than the equivalent 
methane flames (with (H)/(C) =4). An explanation for the 
shrinkage in the observed steady-state zero-g flame can be made 
on the basis of the partial oxidation process. As shown in 
Figures 20 and 21, several complete combustion chemical models 
were used in making both normal-g and zero-g predictions for 
a low fuel flow rate (Q = .22 cc/sec) propylene-air flame. The 
theoretical flame lengths using the various models were: 


Chemistry Model 
Complete 
combustion 


L f (l-g) (cm) 


.60 


L (0 -g ] 


.83 


CO and H 2 

H^O and C(s) .26 


.44 


.35 


Thus, while for a given chemical model the ratio of the normal 
gravity to zero gravity theoretical flame lengths is roughly 
4/3, the ratio can be as low as .58 if one assumes "total 
pyrolysis" in the zero gravity flame, and full complete combus- 
tion for the normal gravity flame. The centerline and radial 
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temperature profiles shown in Figures 22 and 23, and the iso- 
therm map oh Figure 24 show clearly that even using the low 
level energy release CO-H 2 partial oxidation model the tempera- 
ture levels are sufficiently high to promote the formation of 
soot in the fuel rich regions. This conclusion is substantiated 
by the data of Reference 14, which shows that soot can form at 
temperatures of the order of 1000°K depending upon the residence 
time and equivalence ratio. 

Although general qualitative agreement is obtained in terms 
of predicting the increase of flame length and width as the 
gravity level is reduced, it appears that the accuracy of the 
predictions fob the 0-g flame is dependent upon the kinetics of 
the combustion process. 


b. Elliptic Effects 

Although the studies presented so far have indicated a num- 
ber of mechanisms which need additional consideration, their 
importance has been discussed within the context of a parabolic 
system of conservation equations. The question arises as to 
the potential importance of elliptic effects in terms of the 
data that has been analyzed in this study. The significant 
deviations between theory and experiment occur for the lowest 
Reynolds number cases. This has already been partly resolved in 
terms of reduced temperature levels associated with kinetics 
effects and radiation. Nevertheless, axial diffusion effects 
must also play a role in the flame structure at sufficiently low 
Reynolds numbers. For example, boundary layer theory is known 
to apply when the lateral domain of influence is "small" com- 
pared with the axial extent of interest. Specifically, for 
laminar flow the relationship between the lateral dimensions, 

6, and some characteristic axial length is: 



where. 


(47) 



P 


u 

m m 




m 


This choice for Reynolds number is based upon the Dimensional 
Analysis, Section' and (i m and u m are the density, vis- 

cosity, and velocity, respectively, evaluated at the maximum 
flame width point. 

The "condition" for boundary layer theory to apply is: 




(48) 


Although the condition required by Eq. (48) is somewhat 
ambiguous it is sufficient to define boundaries of applicability 
based upon a reasonable requirement for accuracy of the solution- 
What is needed then is a condition that satisfies the require- 
ment that lateral diffusion dominates over axial diffusion- 
For this purpose it is sufficient to require. 


- SI .1 (49) 

L 

for the solution to be "very good." This condition, via 
Eq. (46) implies that the "running" Reynolds number, R e;L * be 
such that: 

R :> 100 (50) 

e L 

Now, since no distinct boundary of applicability can be defined 
a broader definition for the adequacy of the boundary layer 
theory is appropriate. This has been selected according to the 
following criteria: 

R. £ 100 very good 

e L 

25 £ R £ 100 adequate 

e L 

R e £25 poor (51) 
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( 52 ) 


Now, relating R to R e yields the desired limits, viz: 


: L 

R = R 


L 

e r 
m m 


where L is the flame length. 

Figure 25 shows a composite of experimental and theoretical 
results including the boundaries defined by Eq. (51) . The agree- 
ment for normal-g (open symbols) is very good until the “poor" 
domain is encroached. Even there, as already observed, the 
results are actually good, indicating that for normal gravity 
the effect of buoyancy tends to make the basic boundary layer 
criteria conservative. The zero-g comparison is worse as pre- 
viously cited, but the trend is consistent with the boundaries 
defined by Eq. (51). Moreover, the prior analysis of the effect 
of reduced temperature levels including the effects of kinetics 
would result in a leftward shift of the theoretical 0-g curve. 
Thus, better agreement can be anticipated with boundary layer 
theory except at the lowest R e 1 s i.e. - Re m ~ 5, when appro- 

priate kinetics effects and racfiation losses are taken into 
account . 


c. Transient Effects 


Thus far in this discussion emphasis has been based upon 

the importance of kinetics in terms of its relationship to the 

relatively large residence times associated with low Reynolds 

number flows particularly under 0-g conditions. There is yet 

another consideration associated with the residence times 

encountered under zero-g conditions and this is related to 

transient effects associated with drop tower experiments. If 

the non-steady term had been included in our conservation 

equations (Section 3) a parameter characterizing the relative 

importance of transient effects would have appeared. The 

particular parameter is the Strouhal number which relates the 

residence time to the time, T , associated with the event 

c 

which in this case is the drop sequence, viz.: 

r 

Strouhal Number = — ~ (53) 

u T 
m c 
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.Now, referring to Figure 26 it is seen that the character- 
istic residence times, T = r /u are substantially larger 
for the 0-g cases than for the corresponding normal-g condi- 
tions. It is further noted that for T m ~ 100 msec, sustained 
combustion under 0-g conditions is not observed. The data of 
References 1 and 2 show that when quenching occurs it is 
generally preceded by the attainment of a minimum flame length 
in times on the order of 50 to 100 milliseconds. This is an 
unambiguous characteristic time, t g , associated with the 
sudden change in gravity '.level . The Strouhal number, is, there- 
fore, on the order of unity or greater. 

Thus, the hot products of normal-g combustion are retained 
within these 0-g flames because of the extreme reduction in 
convective transport associated with the low Reynolds number 
cases. The retention of the hot combustion gases shields the 
flame zone from fresh oxygen while promoting pyrolysis of the 
fresh fuel being steadily injected directly into the flame 
zone. Thus, before the flame can adjust to a steady state con- 
figuration with convective transport matching diffusive 
transport^ critical reduction in temperature occurs. This 
affects the kinetics of the oxidation process and there results 
a reduced heat release rate until sustained combustion may 
no longer be possible and the flame "quenches." This is indeed 
what is shown by Figure 26. 

It appears then that in addition to the various partial 
oxidation processes described above that for low Reynolds 
numbers, transient effects must also be considered. 

The results of the study thus far have shown that while 
there are mechanisms that need additional attention, the basic 
behavior is predicted for the steady state flames including the 
0-g configurations. Thus, a dimensional analysis has been 
performed and is discussed in the next section. 
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6. DIMENSIONAL ANALYSIS 


A dimensional analysis can serve to (1) define the 
controlling parameters necessary for scaling purposes, 

(2) aid in organizing experimental results and (3) indicate 
the relative importance of the mechanisms controlling the 
process. As cited in Section 2 „ however, the most meaning- 
ful set of values of the flow properties which enter into 
the dimensionless groups are not necessarily the basic initial 
and boundary conditions. In fact* order of magnitude variat 
tions In the transport properties and flow variables are 
common in these diffusion flames. Thus* the approach here 
involves (a) finding a local state within the flow that best 
characterizes the flame and (b) relating those state pro- 
perties to the given boundary conditions. .The result is 
a set of. -relations which' permit the evaluation of the 
Various parameters: including Reynolds, number and Grashof 
number given the boundary conditions. 

Upon a review of the various direct calculations it 
became apparent that the flow field structure is most depend- 
ent upon the properties associated with the heated zones 
within the flow. This suggested that the flame front pro- 
perties shoiald be used including a consistent velocity and 
dimension. The maximum flame width point was selected since 
the maximum flame width radius, r^, is not only unambiguous 
but gives a true characteristic measure of the flame size. 

To verify this hypothesis a systematic study was per- 
formed to determine the implications of this model. Figure 
26 shows a result of correlating the maximum flame width 
velocity, u m , with the associated flame width* r m . This 
theoretical result includes normal-g and 0-g configurations and 
the results indicate that there is a good "universal" correla- 
tion. The significance of this result is twofold: (1) it 

provides one relationship between maximum flame width pro- 
perties and the boundary conditions and (2) it is consistent 
with our earlier representation of the data. The first point 
is discussed in detail later in connection with relating the 
dimensionless groups to the boundary conditions,. ! < 

The second point is crucial in terms of the relationship 
between flame length and the volumetric flow. Examination of 
Eq. (4), Section 3, shows that: 
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2 

Q* u*r* ) 


(54) 


L 


f 


Now, letting the starred quantities be the fuel jet values 
gives : 

L f - Q(~ u q r Q 2 ) (55) 

where Q is the volumetric flow rate of fuel. This is precisely 
what is predicted and what is observed experimentally (c.£. 
Figures 5 and 6). Now from Figure 27 it is seen that: 


2 2 
u r = u r 
mm o o 


(56) 


Thus 0 


L ~ u r 
f mm 


(57) 


which is what Eq. (4) Section 3 „ indicates when the starred 
quantities are replaced by the maximum flame width properties . 
It should be noted that "equations" (55) and (57) may also be 
written : 


and 



(58) 


(59) 


respectively. Thus, introducing the maximum flame width 
properties is thus far totally consistent with the existing 
observations . 
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At this point the study was devoted to evaluating the 
adequacy of the resulting dimensionless groups in character- 
izing the overall flame behavior in terms of the effect of 
buoyancy. It was found that such a meaningful representation 
of the behavior of the methane-air flames in terms of the 
maximum flame width dimensionless groups (Gr/Re^ and 1/Re) 
in explaining the relative importance of the buoyancy, 
viscous and inertia forces in the flame did indeed result. 

For example. Figure 28 shows the theoretical predictions for 
the relative importance of buoyancy, over the range of con- 
ditions of the associated experimental studies in References 1 
and 2 ... These results are consistent with those experiments 
as shown back in Figure 7. Namely, the effect of buoyancy is 
relatively important over the entire range covered by the 
experiments. A composite of the theoretical results is 
shown in Figure 29 showing that at constant Reynolds number 
the effect of inertia increases relative to buoyancy as the 
initial fuel jet radius (and Q) is increased. Conversely, at 
the lower Reynolds number, large r^, conditions, the effect 
of buoyancy is most dominant and these flames are subject to 
the more drastic change upon going from normal-g to zero-g 
conditions. Thus, the effects of transient adjustments and 
kinetics are expected to be most important in this regime. 

Since these effects are intimately connected to the residence 
time an attempt was made to relate the relative theoretical 
effect of buoyancy under normal-g to the Observed quenching 
process. The result is shown in Figure 30. The correspondence 
between the increase in residence time with the buoyancy effect 
is shown and a cutoff for stable 0-g flames is indeed obtained. 
This result is consistent with the results previously shown in 
Figure 26 and, in fact, the correlation shown in Figure 30 will 
carry over when presented as a function of the effective 0-g 
residence time, T m , which in terms of the actual magnitude of 
those times lends further credence to the relevance of the use of 
maximum flame width properties. 

To complete the development of these "effective" character- 
istic parameters the final step relating the maximum flame 
width state to the fuel jet conditions was carried out. 

For a specific fuel/oxidizer/diluent combination, the 
stoichiometric flame temperature, density, viscosity, etc., 
are either known or readily calculated. However, the max- 
imum flame radius, (r m ) and its associated velocity (u m ) are 
generally not known prior to performing a detailed experiment 
or performing a reliable calculation. 
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Hence,, one must relate r m and u m to the specified 
boundary conditions which requires two relations. This 
has been done for the methane-air system with the aid of the 
numerical predictions that have been made during the course 
of this study in comparison with the experimental data 
(Private communication from Mr. Thomas H. Cochram, June 23, 
1971) . 

Regarding the existing data, the essential experiment- 
al data given is the flame lengths and maximum widths and 
whether or not a steady-state zero gravity flame existed. 

Now, the work presented thus far indicates the importance of 
chemical kinetics, particularly for the zero-gravity flames. 
Hence, we were led to correlate the fuel conditions by means 
of a ,E fuei time", t , defined as the burner nozzle radius 
divided by the bulk fuel burner exit velocity. As shown in 
Figure 31, a correlation for r m = r m (T 0 ,r Q ) was indeed 
found. Both theoretical and experimental results are shown 
in Figure 31 and it should be noted that the theoretical 
predictions are within the scatter of the experimental data. 

To complete the definition a second relation was required 
relating r m and u m to the boundary conditions which has 
already been described in connection with Figure 27. 

Now, in Figure 31 it is shown that for normal gravity 
flames (g = 1) it is possible to relate both theoretical 
and experimental values of r m /r Q (within the scatter of the 
results) to t . However, this relation is not unique with 
respect to gravity level. As shown in Figure 32 the relation 
between the average of the theoretical predictions of r m /r Q 
and T q shows a systematic influence due to gravity for non- 
zero gravity levels, and then is abruptly different for zero 
gravity flows. This is indicative of the persistence of the 
importance of the buoyancy effect even for gravity levels as 
low as 0.1-g for the laminar methane-air diffusion flames. 

The relation between r m /r Q and t q is of the form r m /r Q = 
a-m log^QT 0 . Figure 33 shows for non-zero gravity, m is - 

essentially constant, and a may be expressed as a = 10. 5e 
Then, r m /r Q = (10.5 e " e603 ^) - 3.7 log 10 T Q . It is interesting 
to note that the experimental data for methane falls on the 
extension of the straight line shown in Figure 31. The express 
ion for r m /r Q for values of g intermediate between 0 and 0.1 
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can be further defined by more computations. However, here, 
as in any dimensional analysis, we are concerned with order 
of magnitude results in the computation of the dimensionless 
groups, while the variation in r m /r Q would be less than a 
factor of 2 even for very low gravity levels. In summary, 
it appears that the maximum flame width properties best 
characterize the steady state flames that have been investi- 
gated here. In addition, relating the results in terms of 
these parameters has aided in pinpointing the relevance of 
buoyancy and residence time upon the observed flame quenching 
process. 

Now, to provide additional insight into the effects of 
both gravity level and the environmental oxygen concentration 
a parametric study was performed. This study utilized 
Fickian diffusion and "complete" combustion. The results of 
this study are discussed in the next section. 


7 . PARAMETRIC STUDY 


A parametric study can serve several purposes: (1) 

Provide direct predictions on flame structure as affected 
by gravity level and environmental oxygen concentration, 

(2) demonstrate the detail available for comparison with 
in-depth flame measurements and (3) indicate the utility of 
the model do guide the design of experiments tailored to 
provide additional detailed data for verification of more 
refined theoretical models. 

Numerical solutions for the parametric study of the 
two cases were made, and are presented in this section. 
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The initial conditions for/the parametric study were: 


Fuel 

r (cm) 
o 

,0 V 

T K , / \ 

o Q(cc/sec) 

T (°K) 
e 

P (atm) 
e 

Methane 

.051 

294 5.15 

294 

1 

Propylene 

.051 

294 .24 

294 

1 

Fuel 

g 

Environment Oxygen 

Level 

R e 0 

Methane 

i — 1 

A™ 

* 

O 

air, 50%, 100% 


195 

Propylene 

o,h,i 

air, 50%, 100% 


34.5 


In Figures 34 through 43 the theoretical flames shapes 
shown are obtained in the computations made for the above 
initial conditions. Several observations may be made from 
these results : 

For the methane flames, when experimental data were avail- 
able (1-g and 0-g methane/air) there was very good agreement 
between the data and the theoretical flame length and width 
predictions. The initial fuel Reynolds number, Re Q , was 195. 

As shown in Figures 37 amid 38, the effective flame Reynolds 

number. Re , based upon local conditions is over an order of 
m 

magnitude less than R 0 . This behavior is similar to that 
observed by Wohl, et.aS. (Reference 4) in their much higher 
speed butane— air flames; namely, that an initially turbulent, 
or transitional, fuel jet actually developed smoothly into a 
laminar flame. For the methane-air flames, the effective 
flame Reynolds number is still high enough so that a relatively 
long, thin flame develops, even in the absence of buoyancy, as 
shown for the zero-gravity flames of Figure 38. However, for 
the propylene flames (Figures 39 through 43) where the fuel 
Reynolds number is almost an order of magnitude less 
(R e = 34.5), the resulting flame Reynolds numbers are of the 
order of unity, or less, and the flames are very squat. 
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A summary of the theoretical shape factors obtained for 
the methane and propylene flames is shown in Figure 44. The 
relatively stronger influence of oxidizer composition versus 
gravity level upon the shape factor is evident, as is the 
squatness of all of the propylene flames. 

A summary of the influence of fuel and environment com- 
position, and gravity level, upon the buoyancy force is shown 
in Figure 45, expressed in terms of a crossplot of the 
representative flame Grashof number versus the initial fuel 
Grashof number. It may be seen that once the flame survives 
the tendency toward recirculation, typified by a negative fuel 
Grashof number, the most important influence on the strength 
of the buoyant force (and hence the length of the flame) is 
the nature of the environment. Hotter burning oxygen flames 
have higher viscosity levels, and hence mix more rapidly yield- 
ing shorter flames. 

The influence of the environmental composition is further 
illustrated in Figures 46 through 63, showing the axial histor- 
ies, from the nozzle exit to the stoichiometric flame closure, 
of the centerline temperatures, velocities, and viscosities 
of the theoretical flames. The very high temperature, levels 
reached for the pure oxygen and 50% oxygen-50% nitrogen environ 
ments are partly the result of chemical dissociation effects 
(i.e., the formation of OH, H, 0, etc.) not being included in 
the "complete combustion" chemistry model. Plots of the 
centerline viscosity were included to underline the dependence 
of the local, laminar mixing process upon the local temperature 
level and the related combustion process. 

For both the methane and propylene evaluation cases, a 
representative set of conditions was chosen (g = h and a 
50% oxygen environment), and radial profiles of temperature, 
velocity, stoichiometric (equivalence) ratio, Reynolds 
number , Grashof number, and viscosity were plotted and are 
presented in Figures 64 through 74. In addition maps of 
constant mass flow (streamlines) and constant temperature 
(isotherms) were constructed from the appropriate radial 
profiles and are shown in Figures 75 through 77. 

The radial temperature profiles in Figures 64 and 65, 
and the associated isotherm maps in Figures 76 and 77, clear- 
ly illustrate that for the 50% oxygen environment the lack 
of a dissociation mechanism in the "complete combustion" 
chemistry model results in calculated flame temperatures 
that are of the order of 1000 K higher than those that 
would cccur for an equilibrium chemistry model accounting 
for dissociation effects. 
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This excessively high temperature - low density portion 
of the flow field, of course, contributes to the excessive 
squatness of the propylene flame. The effect of using 
"partial oxidation" chemistry models to lower the flame 
temperature level would result in an improved representation 
of the flame chemistry, and in the computation of a narrower, 
more realistic theoretical flow field. However, even then 
the globularity of the low speed propylene flame will result 
in a relatively large surface area, thus contributing to 
potential radiant heat transfer losses. 

In Figure 75, the shapes of typical constant mass flow 
streamlines are plotted. Both the volumetric expansion of 
the streamlines that were initially in the fuel jet as well 
as the entrainment and acceleration of environmental gas may 
be observed. Figures 76 and 77 indicate that high tempera- 
tures will persist on the axis for a signficant distance 
downstream of the stoichiometric flame edge closure. 


8. SUMMARY OF RESULTS 


This report represents a theoretical study of the effects 
of gravity and environmental composition upon the structure of 
laminar, hydrocarbon diffusion flames. The model has been used 
to predict experimental observations and to perform parametric 
studies which have produced the following observations: 

1. For complete combustion the model predicts 
increasing flame length and width with 
decreasing gravity level and decreasing 
environmental oxygen concentration. The 
predicted effect of g-level is observed in 
the bulk of the experimental data whereas the 
effect of oxygen concentration needs verifi- 
cation. However, including partial oxidation 
and reduced temperature levels in the model 
yields predictions which show that a reversal 
in the trend with g-level can occur and may 
be an explanation for some isolated observa- 
tions made on the heavier hydrocarbon. 



2. The accuracy of the predictions is generally 
excellent for normal-g flames while for the 
corresponding steady state zero-g flames 

the predictions at the lower Reynolds numbers 
showed generally wider flames than were 
observed . 

In arriving at these overall findings the importance of 
certain mechanisms controlling the detailed structure of 
diffusion flames have been delineated: 

3. Variable transport properties including non- 
unity Schmidt and Prandtl numbers must be 
taken into account. 

4. A dimensional analysis shows that the state 
of the flow at the maximum flame width point 
characterizes the flame behavior. Further- 
more, this local state can be related to the 
boundary conditions . 


5. The dimensionless groups that characterize 
low speed diffusion flames are the Grashof 
number, Gr, Reynolds number. Re, Prandtl 
number, Pr and Schmidt number. Sc. It is 
further noted that combining Gr and Re to 
form the Froude number, (Fr a (Gr/Re^ ) ~1) 
does not eliminate Re as an independent 
controlling parameter. 

6. Fickian diffusion is shown to be adequate 
although preliminary application of an 
extended model including multi— component 
diffusion provides some improvement in the 
agreement between predictions and experi- 
mental data. 
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7 . 


Shifting equilibrium in the limit of complete 
combustion is an adequate model for the bulk 
of the norma 1-g flames but it is shown that 
partial oxidation models provide an improve- 
ment in the accuracy of the predictions under 
zero gravity conditions. This study suggests 
that partial combustion associated with the 
oxidation kinetics is important. Pyrolysis 
and soot formation are also shown to be of 
potential relevance in terms of their 
effect upon the local molecular weight 
(buoyant force) and radiation. 

8. Incipient recirculation is predicted for 
heavy hydrocarbons when sufficiently large 
negative buoyancy arises. Unfortunately, 
the available data are not sufficiently de- 
tailed to verify that this occurs. However, 
the effect of pyrolysis cited above would tend 
to reduce the magnitude of the negative 
buoyancy and minimize the occurrence of recir- 
culation . 

9. Axial diffusion of mass, energy and momentum 

(elliptic effects) becomes important at very 
low Re, zero-g conditions. This is evident 
from the observation of globular flames, 
particularly as quenching is approached. 

In general, the model presented here has provided a tool 
for accurate predictions of diffusion flame structure over wide 
ranges of conditions including gravity level and environmental 
composition. Furthermore, an improved understanding of the basic 
mechanisms has been gained and regimes of operating conditions 
requiring the inclusion of additional effects have been delin- 
eated. 
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9. 


APPENDICES 


APPENDIX A - FINITE DIFFERENCE SOLUTION 
OF THE CONSERVATION EQUATIONS 


Summary 


The governing flow field conservation equations are 
transformed, into the von Mises plane, and then solved 
numerically employing an explicit finite difference technique. 

Calculation Procedure 

The solution of the system of conservation equations 
presented in Section 2 provides the details of the flow field 
including the velocity, temperature, and species fields. 

Gross characteristics including combustion lengths and flow 
deflections are also obtained. 


The global continuity equation can be eliminated from the 
system of differential equations by introducing the von Mises 
coordinates as the independent variables. The transformation, 
x, r , ~>x, >£, is defined according to the relations: 


N N 
our ~ SI> 'F 


(A-la) 


where 


N N 

- pvr = ^ $ 


x 


(A-lb) 


0 - plane two-dimensional flow 

N = / 

(l - axisymmetric flow 


Introduction of (B-la) and (B-lb) into the differential equa- 
tions results in: 


Element Conservation 

d a i 1 B_r Le gpu 2N ^ a i ~; 

Bx ^N L p r r (A-2) 
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Momentum 


Bu 1_ B P pLou 


2N 


Bu 
B 


i] 


i 

pu 


(p e -p)s 


(A-3) 


Energy 


BH _ 1 B r ^pu 2 n 1 r BH 

Bx ^ Bty L ^ r Pr l B^ 


+Sh, (Le-1) 
i 



g 


(A-4) 


The physical, r, coordinate is obtained by the inverse trans- 
formation : 


N+l 

r 


(N+l) 




ay 


(A-5) 


and the transverse component of velocity, v, is given by: 


v = 




x 



( A-6 ) 


Boundary Conditions 

The governing equations are parabolic and require initial 
conditions at x = 0 and boundary conditions at = « and ^ = 0. 
The initial and boundary conditions are: 


at 


x = 0 , 


0 £ ^ 

3 


u 


u 


o 


H 






x = o ; ^ a . = (a . ) 

J j j e 

u = u (0) , H = H 
e e 


(A— 7 ) 
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rsJ /v 


at x ^ 0 ? $ — » 00 7 a . = (a . ) 

D D e 


u = u (x) , H = H 
e e 


where 


dp 

ST " *>e g = 0 


and 


H = constant 
e 


(A -8) 


(A-9) 


The conditions expressed by Eqs . (B-7), (B-8) and (B-9) with 

symmetry at = 0 completes the specification of the 
boundary conditions. 


The solution of the governing system has been obtained 
employing an explicit finite difference technique (B-l).* 
Figure (B-l) shows a generic point, (n+l,M) in the x->T grid 
network. The finite difference formulation for the calcula- 
tion of the flow at the point (n+l,M) is obtained by using 
the following explicit difference relations where P is any one 
of the three pertinent variables, u, a, or H: 


P - P 

IP. _ n+1 ,M n , M 

dx Ax 


CA-10) 


clP _ 1_ P n,M+l P n,M-l 
2 A'l 


(A-ll ) 


5 r sp "1 

sst L 


b 


n ,M+ 


i, [p 
"2 L n , 


M+l 


- P 


n ,Mj 


b T 

- n , 


n,M-l~* 


A^ 


(A-12) 


where 


b 


2N 

nur g 



(A -13) 


* (B-l) : Zeiberg, S., and Bleich, G., "Finite Difference Cal- 
culations of Wakes," AIAA Preprint 63-449. Also GASL TR-338, 
February 1963. 
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(A-14) 


and 


b n , M+h = h 


f b 

L n. 


+ fc> 

M n , M+l 


] 


* = M(A¥) 


(A-15) 


The conservation equations in difference form are 


Elements 


M = 0: 


l-N. 


1 =( ~ S . 2 (1-N) ax riflu] m i 

a j n+l,o a j n , o (A® 8 L Pr J n,o 


[ (S j ) n.l - ( “j ) n,o] 


(A-16a) 


M^O: 


(a.) ,, = (a.)_ „ + 


M. 


1 Pr 


j n+ l "I'*'* m N (A® 2+N 


n,M.+h (a j >n ' M+ l 


“1 


f* / he / Le b . \ ~ \ • 

L Pr n,M+% + Pr n,M-y a j n,M 


Leb \ 1 

+ Pr n,M-^ a j n,M-l J 


(A-16b) 


Momentum ; 

M=0 


U 


n+l , o 


U + 
n , o 


2 (1+N) Ax 

(a'K 8 


r 1-N "1 r "i 

(pu) a u -u j+(p -o )g 
L p Jn/oL n , 1 n,oJ M e o 


n, o 


(A-17a) 
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M^O 


U 


U + 


Ax 


{(b) i U 

i ' ' ri M4-k T 


n+1 ,M n,M m N ( A ^) 2+n1 - n ^ M+ ^ n,M+l 


fb , + b ,1 U + b Afll U A _)- 

L n,M+^ n,M-h n,M n / M- J 2 n,M-lJ 


+ (o ~p) g — — 
e (nu) 


(A“lVb) 


' n , m 


Energy 
M=0 s 


H 


n+1 


\.o * I ‘.u ) 1 -*!. k >„.= [B n.l- H n.o3 + 


* E< \ i fr _, n.o [fa i , n .l' '“i’n.o 3 } 


(A-18a) 
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H H . + 


Ax r ,b 


n+1 ,M n.M M N (A ^) 2+N 


( (— ) 

\ Pr n. 


H 

M+^ n,M+l 


-[ 


(— ) ! + 

Pr n,lA+h 


,b 


+ 3 Am + n ,M-^ H n,M-l 


+ £ 1 bh . j 

i L 1 pr A, 


(a . ) 


M+h 


i n ,M+1 


? [' bh i + (bh i TTVm-% 1 (a i ) n.M + 


Pr n,M+^ l Pr n,M- 




, ^ n, be ” 1 

+ £ bh. 


, L r Pr •; 
i u “n,M-i§ 


L n,M-l 1 
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The step size in the explicit finite difference scheme is 
controlled by a stability criterion and from studies of linear 
parabolic partial differential equations there results the 
following condition. Reference (B-2): 


r - ] 

(1+N) 6 LLe J n , o S Ax 


£ 


1 

3 



,Le, Le 

# — + ( _Jb ) 

Pr n,M+% Pr n.M-^ 


,(A-19) 

Although the partial differential equations are non-linear, the 
present explicit difference formulation results in a locally linear 

system and Eq. (B-19) provides an estimate of the stable step 
size. The computer program has as an input an arbitrary fraction 
which can be chosen to cut the above step size in the event a 
stability problem arises. 



Figure A-l. - Schematic of the Grid Network Used in the Explicit 
Finite Difference Technique. The Arrows Indicate the Calcula- 
tion of the Flow Field Point n+l,M from Data at Station n. 

(B-2) Richtmyer, R. D., Difference Methods for Initial-Value 
Problems , Interscience Publishers, New York, 1957. 
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The finite difference form for the multicomponent diffus- 
ion analysis , presented in Section 2 , is formulated in a 
similar manner, with the significant addition that the 
diffusional mass fluxes, J^, in Eq. (19 ) must be determined by 
the following matrix techniques; 


- — 


— 

r 

“! 


Sa, 

E „ E, „ E, 


iT 

’ ! F n _j F, 0 . F n . 


1 

1,1 1,2 1 ,N 


1 

| 1,1 1,2 1 , N 

i • 1 

i ■ i 


Sr 

E„ _ E „ „ „ „ s E_ 



hr] 

*1 
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• 

ft 

ft 

*1 
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Sa 2 

2,1 2,2 2 ,N 


2 

' 

i 

i 

2,1 2,2 2 ,N 


3r 

i 

E E . o > E 


J „ 


F „ F n . . . F „ 

! ' act N~i 

N-1,1 N-1,2 N-l/N 


N-l 


N-1,1 N-12 N-XN 

\ Sr 

f 

. , 

P P o o © P 


J N 


0 0 . . . . 0 







1 Sr 

_ — 


_ 


— — 

h- - 1 


Since the solution of the conservation equations is by an 
explicit technique, all of the following parameters which are 
required to solve for the J.'s are known; 


a. - mass fraction of species 
1 

W. - molecular weights 

i 

p - density 


Sa . /Sr 

i 


- the partial derivatives of species in the radial 
direction 


A . 
1 3 


- the diffusion coefficients which form an (NxN) 
symmetric matrix, i.e.. 



The solution vector J comes from 
[E] . [J] = [F] 


(A— 2 0 ) 


48 



The method of solution is by elimination using the largest 
pivotal divisor. At each stage of forward elimination, an 
interchange of rows is performed if and when necessary, to 
insure division by a sufficiently large non-zero element. The 
forward solution obtains the variable in N forward stages. 
The backward solution for each of the remaining variables, 

J N-1" J N-2 ' » « - J j|_ i s obtained by successive substitutions. 
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APPENDIX B-r COMBUSTION MODEL 


The basic model formulation for the laminar diffusion 
flame requires an accurate definition of the thermal field 
rather than the details of the chemical species. Thus, the 
selection of a chemical system in terms of the number and 
kinds of species is somewhat arbitrary. 

Nonetheless equilibrium as well as kinetic computations 
require the specification of species. For a system in 

chemical equilibrium the two principal parameters controlling 
the accuracy of the assumed set of species are the initial 
reactant temperature and the fuel/air ratio. Of course, for 
the determination of the equilibrium state the pressure must 
also be specified. However, the pressure effect is of 
secondary importance. Basically then, the temperature and 
equivalence ratio control the equilibrium state in terms of 
the degree of dissociation. 

In a hydrocarbon-air system the principal products of 
combustion are HO and CO . Furthermore, for fuel lean condi- 
tions 0 will appear while in the fuel rich regime combustibles 
will be 2 present. However, dissociation of the major species 
results in the presence of additional species including H, 0, 

OH and CO, in the fuel lean regime, and under fuel rich 
conditions even more species, including a variety of hydro- 
carbon fragments, C x Hy, and soot, C(s) will appear. At 
temperatures above about 2500°K (which is beyond the experi- 
mental flame temperature.; in air environments ) the dissociation 
of nitrogen will start to become significant and species like 
NO will appear. 

As cited above, the principal quantity of interest is 
the temperature. The objective then was to establish the feas- 
ibility of effectively representing the chemical system for 
the prediction of temperature with less concern for the minute 
details of the species . It was shown in Section 2 that this 
indeed could be done for the conditions of interest in the 
present work (T < 2500°K) . In particular, the entire spectrum 
of -equivalence ratios is modeled by three distinct regimes 
defined according to the fuel oxygen ratio. Figure B.l. 
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In the lean regime complete oxidation of the fuel forming 
CO^ and H^O is assumed. The upper limit for this regime is the 
stoichiometric point where the atom balance is given by: 





+ 2 $ 
'c 



(B.l) 


where 0* s are the molar concentrations. In terms of element 
mass fractions this relationship is given by: 
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2 W. 
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(B.2) 


Accordingly, the lean regime is defined within the limits: 


V i 

W 0 S 2 





^ 0 


(B.3) 


and the specie mass fractions are given by 
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1 a H2 ° 
a H_0 2 H W 
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----- n 
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(B . 4) 


(B.5) 


(B.6 ) 


In the fuel rich side it is observed that CO appears in 
substantial quantities, and depending upon the degree of rich- 
ness, pure fuel appears. 
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In zone A of Figure B-l,it is found that the molar 
concentration of water is essentially constant and that 
carbon is oxidized to CO^, an< 3 CO. This domain is bounded 
by the limits; 


2a . a TT a a , a__ 

c 1_ H _o __c + 1_ H 

W + 2 W rT > W ^ W 2 W TT 
c Hoc H 


(B * 7) 


wherein the species mass fractions are given by 
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In zone B pure fuel, C H , begins to appear and the mole 
fraction of CO decreases. 9ence the bounds are; 
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and the species mass fractions are given by 

a *CH 
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(B . 14) 
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Also note that since nitrogen is assumed not to react in any 
one of the three regimes it satisfies the following relation- 
ship throughout: 




(B . 15 ) 


As previously cited this model represents the thermal 
quite well for the temperatures of interest (T<2500°K) . 


field 
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APPENDIX C 


LIST OF SYMBOLS 


Nomenclature 


P 

Gr 


g 

PI 

h 


h . 

1 


coefficient in the finite difference formulation 
defined by Eq. 

constant pressure specific heat 
binary diffusion coefficient 
Grashof number 
acceleration due to gravity 
mixture stagnation enthalpy 
mixture static enthalpy 

th 

static enthalpy of the i species 
element diffusional mass flux 


1 

k 

L f 

N . 

P 

Pr 


specie diffusional mass flux 

thermal conductivity 
flame length 

number of species 

pressure 

Prandtl number 


Q 


flow rate 


R 

Re 

r 


universal gas constant 
Reynolds number 
radial coordinate 
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half radius 



S flame shape factor 

Sc Schmidt number 

T temperature 

u axial component of velocity 


v 

W 


w . 

i 


W. 

1 


x 


radial component of velocity 
molecular weight 

th 

chemical production of the i species 

th 

molecular weight of the j species 

til 

molecular weight of the k species 
streamwise coordinate 


Q£ 

SKJ 

Gt 


0 

0 

Q 



mass fraction 
element mass fraction 
viscosity 
density 

stream function defined by Eq. (A-l) 
reduced collision integrals 

weighted inverse square of the collision cross section 
collision cross section 
potential energy function 
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kj 

r k j 

^kj 


n 


kj 


x k 


m 


6 

<P 

cr. 


a j 

8 . 

i 


reduced temperature 
molecular weight function 
molecular weight ratio function 
weighted conductivity function 
molecular weight function 
partial conductivity function 
partial conductivity function 
characteristic test time 
characteristic residence time 

boundary layer thickness 

n . . , (fuel/air) 

equivalence ratio ( fuel /air ) 


-) 


stoichiometric' 


th 


collision diameter of k specie 

collision diameter of j *'* 1 specie 

th 

mole fraction of the i specie 


Subscript 
e . 

£ 

f 

i # k 


j 


m 


environment 
center line 
flame 

th 

i and k species 
j*** 1 : element . 
maximum flame width 
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(n,m) generic point in the finite difference grid 

o initial fuel conditions 

t turbulent 

* characteristic value 


specie 

element 


Superscripts 
th 


k 

l 


th 
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Figure 1. - Schematic of vertical laminar jet. 
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Figure 2. - Chemical combustion times for methane-air 
mixtures at 1 atm. 
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Figure 3. - Comparison of quasi-domplete combustion model with 
equilibrium chemistry for a propane-air system. 
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Figure 4. - Comparison of quas i-complete combustion models with 
equilibrium chemistry for a propane-air system. 
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Figure 9. - Comparison of theoretical model with experimental zero-gravity 
flame width in ambient methane-air flames as a function of fuel Reynolds 
number. 
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Figure 10. - Propylene/Air - Max. Q - normal gravity flame 
shapes. 




Figure 11. - Normal gravity tnethane/air flame axial 
cp distribution. 


70 



Oxygen Equivalence Ratio (<p) 



Figure 12. - Zero gravity methane/air flame axial <p distribution 
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Figure 14. - Theoretical flame shape prediction for an 
ambient turbulent hydrogen-air normal— g flame. 
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Figure 15. - Theoretical flame shape prediction for an ambient 
turbulent propane-air normal-g flame. 
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Figure 17. - Comparison of theoretical model with experimental zero-gravity flame 
widths in ambient methane-air flames as a function of fuel Reynolds number. 




Figure 18- - Effect of chemistry model upon analyticsl 
prediction of an extinguished zero gravity flame. 
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Figure 19. - Effect of chemistry model upon analytica 
•prediction of a steady state zero gravity flame. 
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oo Figure 22. - Propylene/Air - min. Q - zero gravity centerline 

temperature histories. 
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Figure 23. - Propylene/Air - min. Q - zero gravity radial 
temperature profiles for the CO partial oxidation 

chemistry model. 
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Figure 24. - Propylene/air - min. Q - zero gravity CO partial oxidation 
chemistry model isotherm map. 
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Figure 25. - Comparison of experimental and theoretical flame 
shape factors as a function of the characteristic flame 
Reynolds number for normal and zero gravity methane-air flames 
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Fiqure 27. - u /u as a function of r /r _ for laminar methane- 
~ mo mo 

air-flames. 
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Figure 28. - Buoyancy -inertia parameter versus viscous-inertia 
parameter as evaluated at the maximum radius flame edge in 
the theoretical normal gravity methane-air flames. 
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Figure 29. - Ratio of buoyancy and viscous parameters 
evaluated at maximum flame radius versus initial fuel 
radius . 
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Figure 30. - Stability criterion for 0-g methane-air flames as a function of 
(G /R ) and r for the GASL theoretical model. 
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Figure 31. - r /r as a function of r in normal gravity methane-air flames. 
■ mo o 
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Figure 33. - The influence of gravity level upon the 
theoretical relation between r /r and r . 
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Figure 34. - Methane-air evaluation case-effect of gravity 
upon theoretical flame shape. 
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Figure 35. — Me thane — 50% oxygen evaluation case effect of 
gravity upon theoretical flame shape. 
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Figure 36 . - Methane -oxygen evaluation case-effect of 
gravity upon theoretical flame shape • 
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Figure 37. - Normal gravity methane flame evaluation 
effect of environment upon theoretical flame shape. 
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Figure 38. - Zero gravity methane flame evaluation cases 
effect of environment upon theoretical flame shape. 
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Figure 39 . - Propylene-air evaluation case-effect of gravity 
upon theoretical flame shape. 












Figure 42. - Normal gravity propylene flame evaluation case 
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Figure 47. - Centerline velocity as 
an air environment. 
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Figure 48. - Centerline viscosity as a function of gravity level for methane in 
P an air environment. 
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Figure 49 . - Centerline temperature as a function of 

gravity level for methane in a 50% oxygen -50% nitrogen 
environment. 
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Figure 50. -Centerline velocity as a function of gravity 
level for methane in a 50% oxygen-50% nitrogen 
environment. 
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Figure 51. - Centerline viscosity as a function of gravity 
level for methane in a 50% oxygen-50% nitrogen environment 
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Figure 52. - Centerline temperature as a function of gravity 
level for methane in an oxygen environment. 
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Figure 53. - Centerline velocity as a function of gravity 
level for methane in an oxygen environment. 
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Figure 54. - Centerline viscosity as a function of gravity 
level for methane in an oxygen environment. 
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Figure 55. — Centerline temperature as a function of gravity level for propylene in an 
air environment. 
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Figure 58. - Centerline temperature as a function of gravity 
level for propylene in a 50% oxygen-50% nitrogen environment 
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Figure 59. - Centerline velocity as a function of gravity 

level for propylene in a 50% oxygen-50% nitrogen environment 
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Figure 60. - Centerline viscosity as a function of gravity 
level for propylene in a 50% oxygen-50% nitrogen environme 
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Figure 62. - Centerline velocity as a function of gravity 
level for propylene in an oxygen environment. 
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Figure 66. - Methane. - 50% oxygen evaluation case velocity 
radial profiles. 
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— Propylene —50% oxygen flame radial velocity profiles 


Figure 67. 




Figure 68. - Methane -50% oxygen evaluation case stoichiometric 
ratio radial profiles. 
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Figure 70* - Methane-50% oxygen evaluation case Reynolds 
number radial profiles. 
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Figure 71. - Propylene -50% oxygen flame radial local 
Reynolds number profiles. 
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Figure 72. - Methane-50% oxygen evaluation case Grashof 
number radial profiles. 
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Figure 73. - Propylene-50% oxygen flame Grashof number radial 
profiles. 
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Figure 74. - Propylene -50% oxygen flame radial viscosity profile 





Figure 75. - Methane - 50% oxygen evaluation case flow 
field streamline map. 
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Figure 76. - Methane - 50% oxygen evaluation case flow field isotherm map 
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